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ABSTRACT 

Relativistic stars are endowed with intense electromagnetic fields but are also subject to oscil- 
lations of various types. We here investigate the impact that oscillations have on the electric 
and magnetic fields external to a relativistic star in vacuum. In particular, modelling the star 
as a relativistic polytrope with infinite conductivity, we consider the solution of the general 
relativistic Maxwell equations both in the vicinity of the stellar surface and far from it, once 
a perturbative velocity field is specified. In this first paper we present general analytic expres- 
sions that are not specialized to a particular magnetic field topology or velocity field. However, 
as a validating example and an astrophysically important application, we consider a dipolar 
magnetic field and the velocity field corresponding to the rotation of the misaligned dipole. 
Besides providing analytic expressions for the electromagnetic fields produced by this con- 
figuration, we calculate, for the first time, the general relativistic energy loss through dipolar 
electromagnetic radiation. We find that the widely used Newtonian expression under-estimates 
this loss by a factor of 2-6 depending on the stellar compactness. This correction could have 
important consequences in the study of the spin evolution of pulsars. 
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1 INTRODUCTION 

The study of the electrodynamics of relativistic stars when these are undergoing normal-mode oscillations has the prospect of providing 
several pieces of information which are relevant not only for the astrophysics but also for the physics of these objects. Astronomical obser- 
vations indicate, in fact, that compact relativistic stars are endowed with extremely intense magnetic fields which reach surface strengths of 
the order of < 10 10 G in older neutron stars associated with recycled pulsars and low-mass X-ray binaries. Observations of young neutron 
stars, on the other hand, show surface magnetic fields that are much stronger and usually of the order of 10 11 — 10 13 G. In addition to this, 
the phenomenology associated with soft gamma repeaters and anomalous X-ray pulsars suggests that the surface magnetic fields can become 
even stronger and up to 10 1 4 — 10 15 G (M ereghetti & Stella 1995, Kouvelio tou et al. 1998, 1999), where they could be produced through 
efficient dynamo processes fThompson & Duncan 1995; Bonanno et al. 2003). 

However, intense electromagnetic fields are not the only important feature of compact objects. The same astronomical observations that 
indicate the presence of such fields, in fact, also make it manifest that these are associated with objects that are very compact and hence 
with strong gravitational fields. As a result, an accurate description of the electrodynamics of compact objects can only be made within 
a framework in which the relativistic corrections are properly accounted for. The literature investigating electromagnetic "test-fields" in 
curved spacetimes (i.e. fields whose energy-density is sma ll when compared with the average rest-mass energy de nsity) is rather extensive 
and has con sidered both spherically symmetric spacetimes iGinzburg & Ozernoy 1 19641: [ Anderson & Cohenll97 0) and also slowly ro tating 
spacetimes (Muslimo v & Tsvgan|[l992llMuslimov & Hardin Jl997l:lKonno & Koiimall200d) . More recently. iRezzolla et alj JJOOI albh have 
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studied the interior and vacuum exterior electromagnetic fields produced by a rotating dipole moment comoving with the star and started a 
systematic approach to provide analytic expressions for the regions of the spacetime near the stella r surfa ce (i.e. the so-called "near-zone"). 
These results have also been extended beyond the low-frequency approximation by Koii ma et alJ 12003I) . who have resorted to a numerical 
approach to find solutions to the Maxwell equations that are approximate, but valid also at large distances (i.e. the so-called "wave-zone"). 

If it is natural to expect strong electromagnetic fields in the vicinity of relativistic stars, it is equally natural to expect that these stars 
too, just like ordinary stars, will be subject to a variety of oscillation modes that become manifest through the emission of both gravitational 
waves and electromagnetic waves. Because both of these observational windows may soon provide data that could reach us simultaneously, 
they offer unique opportunities to test the physics and internal structure of these objects as well as the properties of matter at nuclear densities. 
The investigation of the electromagnetic fields produced by an oscillating magnetized star is, however, a rather complex problem and this 
is testified by the scarce literature produced so far in this context. Exceptions in this sense are, to the best of our knowledge, the works of 
McDermott et al. (1984) and Muslimov & Tsygan (1986). While the first authors (McDermott et al. 1984) have modelled the electric fields 
produced in the wave-zone by a neutron star subject to both toroidal and spheroidal modes to calculate the electromagnetic damping of these 
oscillations, Muslimov & Tsygan (1986) have instead computed the exact analytical solutions for the electromagnetic fields produced, both 
in the near-zone and in the wave-zone, by an oscillating neutron star with a dipolar surface magnetic field. In both cases the neutron star was 
considered in vacuum and, for simplicity, the stellar magnetic field was considered unperturbed. More recently, Timokhin et al. (2000), have 
dropped the assumption of a vacuum region outside the stellar surface and performed an extensive investigation of the electrodynamics of 
the magnetosphere when this is perturbed as a result of toroidal oscillations. 

An aspect that these three works share is the use of a Newtonian description of gravity and, as we have argued above, this is expected to 
be a poor approximation in the vicinity of the stellar surface. The aim of the present work is to layout a general formalism for the calculation 
of electromagnetic "test-fields" in the curved spacetime of a spherical relativistic star subject to perturbations. We do this by extending 
the approach developed for the electrodynamics of rotating relativistic stars (Rezzolla et al. 2001a, b) to includ e the pos sibility that the 
conducting crust may have a velocity field as a result of the stellar oscillations. As in our previous investigations JRezzolla et all2001albh . 
our main aim here is to provide analytic expressions in a form which is as simple as possible and can therefore be used also in astrophysical 
applications. In this first paper we will concentrate on the development of the mathematical formalism and in particular on the solution in 
the two distinct regions of the space for which exact solutions can be found, i.e. near the stellar surface where the relativistic corrections are 
most important (near-zone), and at very large distances from it where the spacetime approaches a flat one (wave-zone). In both cases we find 
that two main corrections emerge in the expressions for the relativistic electromagnetic fields and these are associated: i) to the amplification 
of the electromagnetic fields on the stellar surface produced by the local nonzero spacetime curvature; ii) to the gravitational redshift which 
will affect the electromagnetic waves as these propagate in the curved spacetime. 

The expressions for the components of the electromagnetic fields presented here for the two regimes are completely general and 
do not refer to a specific magnetic field topology nor to a precise velocity field. The application of these solutions to the case of those 
spheroidal and toroidal velocity fields most frequently encountered in the study of stellar oscillations will be presented in a companion pa- 
per iRezzolla & A hmedov 2004). Here, however, to validate the expressions derived and to provide a comparison with known results, we 
consider a background magnetic dipole and the simplest of the perturbation velocity fields: the one produced by the uniform rotation of the 
dipole. In th e near-zone, the compa ri son of the derived e xpressions for the electric and magnetic fields can be compared with the solutions 
calculated bv lRezzolla et al]l2001al) : lKoiima et alJJ2003ti . while, in the wave-zone, the new components of the electromagnetic fields can be 
used to calculate the general relativistic expression for the electromagnetic energy loss through dipolar radiation. This validating test is also 
useful to r eveal that the well-known N ewtonian expression commonly used to estimate the rotational energy loss through dipolar radiation 
in pulsars l Pacini 1968; Gunn & Ostriker 1969; Shapiro & Teukolskv 1983) under-estimates this loss by a factor 6 if the star is very compact 
or by a factor of 2 if the star is less compact. 

The paper is organized as follows: in Section|2|we formulate the basic Maxwell equations in the exterior Schwarzschild spacetime, while 
in Section[3]we discuss the expressions for the electromagnetic fields in the stellar interior and the boundary conditions that constraints them 
at the stellar surface when this is infinitely conducting and subject to velocity perturbations. Sections|4]and|5]are devoted to the derivation of 
the vacuum electromagnetic fields in the vicinity of the source and at large distances from it. There, we provide the exact analytical solutions 
to the general relativistic Maxwell equations for arbitrary background magnetic fields and velocity perturbations. Section|6|discusses how 
the formalism developed in the previous Section can find a simple and yet useful application in the case of a rotating magnetic dipole 
and provides the general relativistic correction to the well-known expression for the energy loss through dipolar electromagnetic radiation. 
Section^contains our conclusions and the prospects for further applications of the formalism presented here. A number of details about the 
calculations carried out in the main part of the paper are finally presented in the Appendices A-D. 

We here use a spacelike signature (— , +, +, +) and a system of units in which c — 1 (unless explicitly shown otherwise for conve- 
nience). Greek indices are taken to run from to 3, Latin indices from 1 to 3, and we adopt the standard convention for the summation over 
repeated indices. Furthermore we will indicate four-vectors with bold symbols (e.g. u) and three-vectors with an arrow (e.g. u). 
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2 MAXWELL EQUATIONS IN A SCHWARZSCHILD SPACETIME 

The investigation of the response of a relativistic and magnetized star upon a generic perturbation would require, in principle, the use of 
the coupled system of the Maxwell and Einstein equations. Such a formulation, however, is overly complicated for most astrophysical 
applications and a few reasonable assumptions can be made to simplify the treatment to a more tractable form. The first of these assumptions 
concerns the corrections to the Einstein field equations introduced through the mass-energy contribution of the electromagnetic fields inside 
and outside the relativistic star. It is not difficult to show that these electromagnetic corrections are proportional to the electromagnetic energy 
density and are rather small in most compact and magnetized stars. Indeed, if (po) is the average rest-mass density of a star of mass M and 
radius R as measured at infinity, these corrections are at most 

B 2 „„ „„_ 7 / B \ 2 flAM Q \ f R x3 
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when compared with the total rest-mass energy density. Hereafter we will assume these corrections to be negligible for the stellar objects 
which we are interested in and thus treat the electromagnetic fields as "test-fields" in a given curved background. Other authors, interested 
also in the structural modifications produced in very highly magnetized stars (Bocquet et al. 1995, Konno et al. 1999, 2000, Oron 2002, 
Ioka & Sasaki 2004), have preferred not to make this approximation and have instead included the magnetic contribution to the spacetime 
curvature. The second simplifying assumption is that we will consider also negligible the corrections produced by a global rotation of the 
spacetime and induced by the rotation of the compact star (these corrections were discussed by Rezzolla et al. 2001a,b). The numerical 
analyses carried out by Geppert et al. (2000) and Zanotti & Rezzolla (2002) show that this is a rather good approximation in most cases of 
astrophysical relevance. 

As a result of these two assumptions we can work in the background spacetime of a nonrotating star, whose line element in a spherical 
coordinate system (t, r, 6,4>) is given by 

ds = goo{r)dt + gu(r)dr + r d6 +r sin 6d<j) . (2) 

The portion of the spacetime external to the star (i.e. for r > R) is simply given by the Schwarzschild solution with —goo = N 2 = 
(1 — 2M/r) and jn = 1/goo- For the portion of the spacetime interior to the star, on the other hand, the metric functions can be specified 
in terms of two potentials A(r) and $(r) so that 

2$(r) 2A(r) ( \ 2lTl(r)\ 

9oo = -e y ' , gu = e >• > = 1 ^ , (3) 



where p(r) is the total energy density and m(r) = 4-7T J r 2 p(r)dr its coordinate volume integral. The precise form of these potentials is 
obtained through the familiar solution of the Einstein equations for relativistic spherical stars, i.e. the TOV equations (Shapiro & Teukolsky 
1983) and the metric functions are then matched continuously to the external Schwarzschild spacetime so that 

-goo(r = R) = N 2 R = 1 - ?M- , and gil (r = R) = -^ . (4) 

Within the external portion of this background spacetime we select a family of static observers with four- velocity components 

(u a ) obs = N-'f 1,0,0,0) , ( Ma ) obs = n(- 1,0,0,0 ) , (5) 



and associate to them orthonormal frames with tetrad four-vectors e^ = (eg,e#, eg, ei) and 1-forms uj^ — (uj 1) ,uj r , w , w ), whose 
components are 

e£ = -i(l,0,0,0j, wl = N(l,0,0,o\ , (6) 

e? = n(o,1,0,o\, u f a = ±(o,l,0,o\, (7) 

e£ = in),0,l,0j, w« = r(o, 0, 1, J , (8) 

^"^im^f 0,0 ' ' 1 ) ' ui=rano(o,0,0,l\, (9) 

and which will become useful when determining the "physical" components of the electromagnetic fields. 

The relations among the electric E a and magnetic B a four- vector fields measured by an observer with four- velocity u a can be expressed 
through the electromagnetic field tensor F a p 

F a p = 2u [a E l3] + rj a01 su' f B s , (10) 

where T\ a m = 2(^0,3 — Tp a ) and r/ap^s is the pseudo-tensorial expression for the Levi-Civita symbol e a ^js- 
The first pair of the general relativistic Maxwell equations 

3lF [al3n] =2(F aflrf + F- fa ,p + F f 3- l , a ) = 0, (11) 
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can then be projected along the tetrad four-vectors carried by the proper observers and expressed in terms of the "physical" electric and 
magnetic three-vectors E and B (Thorne et al. 1986) as 

V • B = , 9,B = -Vx (NE) , (12) 

where V represents the covariant derivative with respect to spatial part of the metric 0. Note that equations 1121 are not valid in a spacetime 
admitting a rotational Killing vector (e.g. such as the one produced by a rotating relativistic star) and the general form they assume in that 
case is given in Appendix IdI 

Once expressed in their "hatted" component-form, equations <12> simply become 

sinOdr (r 2 B fA j + j-d e (sinflB 5 ) + jjfyB* = , (13) 

;m, ' A a t B f = d<,E?-d e Um9E*\ , (14) 



N 
^^jdtB 8 - -d*E f + S mOdr(rNE^, (15) 

(jjj&B* = -ft. (rJVB*) + d 6 E? . (16) 

Indicating now with w the conductor four-velocity and with p e the proper charge density, p e w a will represent the convection current, 
and the second pair of Maxwell equations is then written as 

F^. fj = 47vr=MpeW a +n, (17) 

with the semicolon indicating the covariant derivative with respect to the metric J5J and J being the total electric-charge current. Also these 
equations can be written in terms of the physical electric and magnetic vectors as 

V-E = 4-K Pe , d t E + 4-Kj = -Vx (NB) , (18) 

or, in component form, as 

N sinOdr (r 2 E'~') +rd e (sin8E e ) +rd< l ,E i = 4np e r 2 sin 9 , (19) 



de (sin 0B*) -d^B 6 = ( ^p J dtSf + 4nr sin f , (20) 

d4,B f -sinOdr (rNB^ = (^r^) d t E § + 4-Krsm9J § , (21) 

d r (rNB § ^ - d B f = (^ dtE* + 4-Kr.f . (22) 

Note that the total electric-charge current in equations I17i includes the conduction current j a , which is associated with electrons having 
electrical conductivity a, and that Ohm's law can be written as 

j a = aFcpw 13 . (23) 

So far, we have considered the star to be static and in equilibrium, but we can now introduce a fluid perturbation in terms of a 4-velocity 
Su, so that the fluid velocity is in general w = 8u, where the 8 indicates that this is an Eulerian perturbation. The components of the velocity 
perturbation are then 

6u- = r(lM)=r(l,e- A 6v*,^,-^) , (24) 

V / \ r rsmOJ 

and 

8u a —ri —e ,e 8v r ,r8v ,r sin 05v J , (25) 

where Sv 1 = dx l /dt is the oscillation 3-velocity of the conducting stellar medium and Sv' = u>k8v k are the components of the oscillation 
3-velocity in the orthonormal frame carried by the static observers in the stellar interior 1 . Because we are interested in small velocity 
perturbations for which 8v' <C 1, we can neglect terms 0(8v 2 ) and use the normalization for the four-velocity w a w a — —1 to obtain 

-1/2 

r = 



8v'8v k 

— ffOO 1 + Qih 



e~* ■ (26) 

Note that Su is the only perturbation which we need to consider here. While, in fact, other perturbations (e.g. in pressure and in energy 
density) may be present and influence the stellar structure, these will not affect the induction equation which will depend linearly on 5u only 
[cf. equations <10H . Furthermore, for simplicity we will assume that Su is a given, generic function and is not necessarily the solution of an 

1 We indicate with a tilde the 1-form basis of the static observers inside the star to distinguish this from the corresponding basis given in |S|. 
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eigenvalue problem (this was considered by Messios et al. 2001). Finally, because we treat the electromagnetic fields as "test-fields", we will 
not consider any perturbation to the spacetime metric J2j, and therefore work within the so-called Cowling approximation. 

Note that the introduction of perturbations in the velocity will produce important modifications in the electric and magnetic fields in 
three different spatial regions: i) the stellar "interior", r < R; ii) the "near-zone ", i.e. a vacuum region outside the stellar surface with radial 
extension r > R and Hi) the "wave-zone", i.e. far away from the star, at r 2> R, where the perturbations propagate as electromagnetic 
waves. Each of these modifications has observable consequences and the following Sections are dedicated to discussing two of these regions 
for the general case of arbitrary velocity perturbations and generic magnetic field topology. 



3 ELECTROMAGNETIC FIELDS IN THE STELLAR INTERIOR 

We here concentrate on the modifications in the internal electromagnetic fields produced by velocity perturbations of different types. While 
not directly observable, these electromagnetic fields, and in particular the values they assume on the stellar surface, leave an important imprint 
on the electromagnetic waves produced by the oscillations and that will reach infinity. 

We here start by using Ohm's law I23i and the four- velocity 1241 . to obtain the explicit components of the total current J a as 

f = p e e^Sv f + o\E f + e"* (dv S B* - 6v* B S \] , (27) 

J 8 = p e e-*8v i + <r[EP + e-*(8v*B f -6v*B* , \] , (28) 

J* = p e e _ *fo^ + a \E $ + e~* Uv^B 8 - 5v S B*\~\ . (29) 

A convenient way of simplifying the problem is that of considering the fluid as perfectly conducting, i.e. with a — ► oo (this is the limit 
of ideal magnetohydrodynamics). While idealized, this approximation is a rather good one in the present case, since the Ohmic diffusion 
timescale is several orders of magnitude larger that the typical timescale for the stellar oscillations. Once this assumption is made, the electric 
field in the interior of the star can be easily derived from <27> - 1291 to be 

EL = -e~* (6v § B? nt - 6v*BL) , (30) 

EL = -e~* (5v*BL - 5v f B&) , (31) 



EL = -e~* Sv'BL - Sv'BL) , (32) 



where we have indicated with an index "int" the field components for r < R. The Maxwell equations 1131 — 1161 combined with the expres- 
sions 1301 — I32i will lead to the following set of the magnetic induction equations 

(r sine) d t Bf nt = d^(dv f Bi t -dv $ Bt nt )+d e [sin9(dv f Bi t -dv 8 Bt llt )] , (33) 

(r sin 9) dtBtt = d (Sv § B? nt - Sv $ bL) - sin 9e- A d r [r (5v*bL - Sv'bL) } , (34) 

rd t B? nt = e' A d r [r [Sv^B[ nt - Sv f Bf nt )] - do (Sv § B? nt - 6v*bQ . (35) 

The boundary conditions for the magnetic field across the stellar surface r — R can then be obtained after requiring continuity for the 
normal (i.e. r) component 

Bext|r=il = B r R , (36) 

while leaving the tangential (i.e. 6 and <f>) components free to be discontinuous through surface currents 

BL\r=a = B{+A3rI $ , (37) 

Bfj^n^Bl-Anl 8 , (38) 

where I k are the components of the surface current and B l R = Bl nt \ T= R (the index "ext" has been here used to refer to fields at r > R). 
In a similar way, the boundary conditions for the electric field across the stellar surface can be derived from imposing the continuity of the 
tangential components, while leaving the normal one free to be discontinuous through a surface charge distribution E s . Simple algebra then 



1 '-..«□* X„.<f> D» 
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gives 

sL|r=i? = EL|r=fl + 4^E a = -^- (*£ B* - Sv R B R ) + 4^E S , (39) 

EL \ r=R = EL I r=* = ~ (to* B R - Sv R B R ) , (40) 

*L\r=R = EL\r=R = ~ (Sv R B § R - 5v R B f R ) , (41) 

where 5v R = 5v z \ r =R are the oscillation velocities at the stellar surface, the only relevant ones in this paper. 



4 ELECTROMAGNETIC FIELDS IN THE NEAR-ZONE 

We now consider the form of the electromagnetic fields generated by stellar oscillations in a region of vacuum space exterior and close to the 
stellar surface. Our treatment in this Section will be as general as possible and we will not restrict it to any specific magnetic field topology or 
velocity field. A more detailed discussion of the electromagnetic fields produced by those spheroidal and toroidal velocity fields most often 
discussed in the physics of compact stars will be presented in a companion paper (Rezzolla & Ahmedov 2004). 

Given a characteristic frequency for the stellar oscillations /o, the spatial extension of the near-zone can then be estimated to be 
A ~ c/fo, which is the distance travelled by an electromagnetic wave in a timescale l//o- For a typical neutron star of mass A/ and radius 
R, the fundamental frequency of oscillation is /o ~ y/GM/R? ~ 2 — 3 kHz, so that the near-zone extends from R to ~ 100 R. Since we 
are neglecting terms 0(5v 2 ) in the Maxwell equations, we need not account for the displacement currents in equations I20i — I22i and set 

d t E f = d t E e = dtE* = , (42) 

so that the second pair of Maxwell equations can be written as 

N sin 6dr (r 2 E f ) + rde (sin 6E e ) + rd^E^ = , (43) 

d e (sin 9B^\ - d^B § = , (44) 

9 B f - sin ed r (rNB$\ = , (45) 

d r (NrB § ) - d e B f = . (46) 



The solutions to the Maxwell equations with sources represented by a bounded distribution of currents that vary harmonically with 
time is mostly easily found if the electromagnetic fields are expanded in their multipolar components. An efficient method to obtain this 
multipolar expansions has been proposed by Bouwkamp & Casimir (1954) in Newtonian electrodynamics and uses the general properties 
of solenoidal vector fields. We here follow the same approach but extend it to relativistic electrodynamics. Consider therefore four scalar 
functions S, T, X, and Z, through which the general solutions to the vacuum Maxwell equations <13> - J16t and J43I — I46I can be written in 
the form 

B f = - - X 2n [sing9 9 (smOdgS) + o^S] , (47) 

r A an 

B § = —OedrS + AT -. a d^Z , (48) 

r Nrsmo 

B" = Z^-a d ^rS ~ TF^-zdeZ , (49) 



and 



r sin Nr sin 



1 



E r = — T^-x sin Ode (sin 6d e T) + diT \, (50) 

r 2 an 

E § = -dgd r T + - -. d^X , (51) 

r Nrsmo 

E* = JL-d^drT - -±-d e X . (52) 

r sin Nr 

We will refer to the functions S, X and T, Z as to the "magnetic" and "electric" functions, respectively. Furthermore, to handle analytically 
the angular derivatives, we will assume they can be separated in variables and expand the angular dependence in terms of spherical harmonics 
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Y tm (9,4>) 

oo £ oo £ 

S{t,r,6,4>) = Y. E Si m {t,r)Yi m {0,cf) , 



X(t,r,8,0) = Yl E X im {t,r)Y h 



1=0 m=-t 



1=0 m=-t 



T(t,r,M) = £ E T irn (t,r)Y im (e,<j>) , Z(i,r,M) = X] E Z em (t,r)Y em (9,cb) . 

1=0 m—-e 1=0 m=-l 

Using this ansatz and suppressing the summation symbols over the £, m indices, the general solutions 1471 — 1521 then take the form 

B *- £C£±l) g v 

E> — ^ Semllm , 

B — — drSe m doYi m + — — : — -Z( m dd,Yi m , 
r Nrsmu 

i N 1 

B — : -drSlmdtiYtm T—Zl m 'deYl m , 

rsintf Nr 

for the magnetic field components, and 



E r = 



e(t + i) 



Ti m Yt, 



E = drTlmdgYlm + — — : — -Xl m dA,Yl„ 

r Nrsvad 



E = — : — -drTlmd^Ylm — — — Xt m d$Yi r , 
rsmt) Nr 



(53) 
(54) 

(55) 
(56) 
(57) 

(58) 
(59) 
(60) 



for the electric field components. We next concentrate on how to use the Maxwell equations to recast expressions 1551 — 1601 in a form more 
useful for astrophysical applications. 



4.1 Magnetic Fields 

Substituting the general solutions 1551 — I57i for the magnetic field into the Maxwell equations 1441 — 146^ gives two equations for the unknown 
functions Z( m and St m 



1 



d_ 
eli- 



de (sm9deY im ) + . 2 d%Y ln 



2M\ d 

i- -TOIirn 

r I dr 



e(£ + i) 







Sim — . 



(61) 



(62) 



which have the important property of depending on Z or S only, so that they can be solved separately. The first equation has indeed the trivial 
solution Ztm = since the content of the square brackets in J61I is a function of 9 and <f> only. The solution of J62t . on the other hand, is 
less straightforward. It is convenient to introduce the new variable 

- i r 

and factor out an explicit quadratic radial dependence in the functions Si m which are then redefined as 



(63) 



St 



2 h(t,i 



Using now 1631 and 1641 . equation <62t takes the form 



d_ 
dx 



{(i±|)A [(1 _ )2 ,]} + , ( , + 1) ^o. 



which can be solved in terms of Legendre functions of the second kind Qe(x) (see Rezzolla et al. 2001a for details) so that 



St m (t,r) - 



(1-xf d 



(1 + x 



dQe 
dx 



se m {t) 



(64) 



(65) 



(66) 



M dx 

Note that all of the time dependence in 1661 is contained in the integration constants stm (t) an d these are determined after suitable boundary 
conditions across the stellar surface are imposed. 

We next consider how these expressions vary when a perturbation, caused for instance by a nonzero velocity field 1241 . is introduced. 
Hereafter, we will assume that a background magnetic field is present and which varies with time only on a timescale which is much 
longer than that of the stellar oscillations and can therefore be considered static. Hereafter we will indicate with B the components of 
this zeroth-order (background) stellar magnetic field and with SB 1 (t) the first-order perturbations, which clearly are the ones possessing a 
time-dependence. 
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Within a linear regime in the perturbations, a convenient way to isolate the perturbed part of the magnetic fields is suggested by the 
structure of the solutions <55> - <57> which confines all of the time dependence in the S£ m terms. As a result, it is possible to express the 
magnetic field perturbation in terms of new, time-dependent integration constants 5sg m (t) that simply add to the background ones Si m . The 
new components of the magnetic field generated by a velocity perturbation will therefore have the generic form 



B r = B r n 



8B r (t) = 



B=Bo+5B a (t) = - - 



Bf = B% + 5B' p (t) = 



e(£ 


+ 1) d 


(1 


M 3 dx 




1 + x 






(1-x) 


3 




1 + x 




(1-x 


) :i J 



(1 + x)^ 

dx 



[sim +SSl m (t)]Ytr, 



1/2 



M s dx* y 



1/2 



1 d 

IPdx' 



dx 



dx 



^^ 



(1+X 



dQ e 

dx 



[stm + 3si m (t)]dgYi rl 



[stm + Sse m (t)] — — -d^Ye„ 



(67) 



(68) 



(69) 



The values of the integration constants Sst m (t) can be calculated rather straightforwardly if the oscillation modes are assumed to have 
a simple harmonic time dependence of the type exp(— iuot), where u> is the mode frequency. In this case, in fact, the radial component for 
the induction equation I33i at the surface of the star becomes 

d t (5Bi) = -iuJB% = ^i-^ [dv (SvlBi - SviBi) + d e [sinS (6v* R B 9 R - Kb r )] } , 



(70) 



1 "■ Rsin9 

and the boundary condition for the continuity of the radial magnetic field 1361 can then be used to determine the coefficients 8si m (i) (see 
AppendixlBlfor details). Note that this procedure is different from the one used by Timokhin et al., (2000), in which the integration constants 
are instead obtained through the boundary conditions applied to the continuity of the tangential components of the electric field. Of course, 
the two procedures are equivalent and yield identical results. 

A few remarks are worth making at this point. The first one is about the assumption of infinite conductivity for the stellar material, which 
implies that the magnetic field is adverted with the fluid when this possesses a nonzero velocity. This is expressed by the general relativistic 
"frozen-flux" condition 2 

d t 5B = V x (5vx Bo) , (71) 

and is, indeed, verified by expressions 1671 — I69i . The second one is that these expressions, although interesting for their generality, are not 
particularly useful if we do not restrict our attention to a specific background magnetic field configuration; this will be done in Section l6~Tl 
where we concentrate on the specific form that equations <67t - j69t assume for a background static dipolar magnetic field. Finally, note that 
no reference is made at this point on the type of velocity perturbations introduced as these are effectively incorporated in the expressions 
for the functions 5s£ m (t); in Sections l6.2l and l6~3l we will discuss the form of these functions for the case of a magnetic dipole in uniform 
rotation. 



4.2 Electric Fields 

We next concentrate on the form of the electric field components and consider for this the solution to the Maxwell equations <14M16l and 
I43i . Hereafter we will assume that the star has a zero net electric charge (i.e. f p e ^/-yd 3 x = 0, where 7 is the determinant of the three-metric 
and d 3 x is the coordinate volume element) and that the background magnetic field is stationary (i.e. dtB^ = 0). In this case, the background 
electric fields will be identically zero (i.e. Eg — 0) but time-dependent, first-order electric fields SE l (t) may be present and induced by the 
perturbations in the background magnetic field. 

To derive the expressions for these perturbed electric fields we substitute the decomposed expressions I50i — 1521 into the Maxwell 
equations 1141 — 1161 and 1431 to obtain the following set of equations for the unknown functions Ti m and Xi m 



d_ 
dr 



1 _ ™L) ± Tl , 
r j dr 



e{e + i) 



Tim — . 



(72) 



1(1 + l)Y em Xer, 



-r 2 d t B f 



(73) 



drXlmdeYim — -—d t B 



drXlmd^Ylm = — — dtB 



N 



(74) 
(75) 



2 Note that the lapse function does not appear multiplying the magnetic field in expression 1711 . This is because of a cancellation when Ohm's law is used 
to replace the contribution to the Maxwell equations coming from the electric field (see Appendix |D] for the expression of Ohm's law in terms of physical 
three-vectors). 
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It is apparent that equation 1721 is the same as I62t but for the unknown functions Tt m and will therefore have solutions in terms of 
Legendre functions of the second kind Qt(x) and the integration constants 5t( m {i) 



T tm = {l-x) 2 M 2 ^- 
dx 



C + "T 



5t lm (t) . 



(76) 



In general, that is for a multipolar magnetic field, the "electric" functions Xp m (i) are determined directly from the time derivatives of the 
"magnetic" functions Si m [cf. eqs. J55I and 1731 1. i.e. 



Xtm — —dtStm — 



{l-xf d 






Sxtm{t) 



(77) 



M dx 

where 8xt m (t) — — dt[5se.m{t)]. Expression 1771 underlines that the functions Xi m represent the contribution to the electric field coming 
from the time variation of the magnetic field, with the coefficients dxe m being determined, also in this case, through the use of suitable 
boundary conditions. 

Using expressions 1761 and 1771 and standard recurrence formulae for the Legendre functions that we will recall in Appendix A, the 
solution for the vacuum electric field 1581 — 1601 can be written as 



SE r 



f (I + 1) 
(1-z) 2 



[Ql-! -(i+l-ix) Qt] 5U m {t)Ylr, 



5E e = - y Tn„ [xQi - Qe-i] 5t em (t)d e Y im + v 



dE* = 



(l-x) 2 N 

e (i + 1) 



[xQe - Qe-i] 8tt m {t)d<f,Ytm - 



N sin 9 dx 

l-x)M d 



ft i ^ d ® 1 



'1 + x) 



8xirn{t)d<i>Ye 
dQ ~ 



dx 



Sxe m (t)dgYem , 



(78) 
(79) 
(80) 



(l-a;) 2 iVsin0 L ""^ c ~ s *-~ iJ "*'" ^'~*~"'- jy dx 

where, again, Qt — Qt(x). Note how equations 1781 — 1801 clearly indicate that given a magnetic field with multipolar components up to the 
order t, the corresponding electric field will have multipolar components up to the order £ + 1. This is another aspect of the interconnection 
between electric and magnetic fields contained in the Maxwell equations. 

Finally, the integration constants 8te m and Sxi m can be calculated from the jump conditions 1401 — 1411 of the tangential electric field 
1791 — I80i at the stellar surface and are given by 

(1 — X ) 2 

St tm(t) = ^ rn R ' [x R Ql{x R ) - Ql-l{X R )] X X 



P (£ + 1) 



8x irn (t) = --- 



x J dn \d e YL [svi(t)B a - 5vl(t)B R ] - *^p [*£ (t)s£ - <5<(tK 



dQ d e Y? m [Sv 6 R (t)Bl - 5v R (t)B R ] + ^Jf [Svt(t)B R - Sv R (t) B i 



Mf R j i l - • • -J sine 

where Qi(x R ) = Qi (1 — R/M), dfl = sm9d9d<f>, and f R is just shorthand notation for 



fa 



3 ( R_ 

M 



lnNt + 



2M 



1 + 



M 



3 f B_ 
"8 Im 



In 1 



2M\ 2M , 



M 
~R~ 



(81) 
(82) 

(83) 



5 ELECTROMAGNETIC FIELDS IN THE WAVE-ZONE 

The treatment of the electromagnetic fields in the wave-zone can start from considering the Maxwell equations 1131 — 1161 and <19> - <22t in 
the case when the electric charges and currents are absent. In particular, inserting equations I15i and <16> in the time derivative of equation 
J20I and using the expression for 1191 in vacuum (i.e. for p e — 0), one obtains the following wave equation for the radial component of 
electric field on a curved Schwarzschild background 



dE r 



N 



>[ 



N 2 d r r 2 E 



N' 2 



d e sin 9d e E 1 






(84) 



r A sm a \ / r" sm* 

Substituting the decomposition J58> for the radial component of the electric field in terms of spherical harmonics, into the wave equa- 
tion I84i yields the well-known Regge- Wheeler equation describing the dynamics of vector perturbations in a Schwarzschild spacetime 
(Chandrasekhar, 1992) 



N 2 



(dt - d 2 r , ) T lm + lit + 1)—T t 



0. 



(85) 



where r* = r + 2M In (r/2M — 1) is the tortoise coordinate (Similar equations can be obtained also for the other functions Se m ,Xt m , and 

Zlm-). 

Although equation <85l has been extensively studied in the past and contains information about general relativistic effects such as the 
scattering off the curved background spacetime (see Malec et al. 1998 for a recent investigation), it does not represent the most convenient 



10 L. Rezzolla, B. J. Ahmedov 

way to evaluate the components of the magnetic and electric fields in the wave-zone. There, in fact, the relativistic corrections 0(M/r) can 
be reasonably neglected and we can seek guidance in determining the components of the electromagnetic fields from simpler Newtonian 
expressions. 

We note that assuming a flat background spacetime for the form of the electromagnetic fields in the wave-zone is clearly an approxi- 
mation, but it is not in contrast with the general relativistic approach followed so far. As will become clearer in the following, in fact, the 
final expressions which we derive on a flat background, still need to be completed through the specification of suitable boundary conditions 
to be imposed at the stellar surface and these will indeed contain the general relativistic character of the curved spacetime. Furthermore, 
adopting a Newtonian framework as a guideline is not only a good approximation, but also has two important advantages. The first one 
is that is possible to find analytic expression to the electromagnetic fields which will depend on outgoing spherical Hankel functions and 
their derivatives [cf. equations J97I — I102l l. In a fully general relativistic treatment, instead, the numerical solution of an ordinary differential 
equation is necessary to obtain the form of the electromagnetic fields in the wave-zone (this will be discussed in Rezzolla & Ahmedov 2004). 
The second advantage is that in the specific but important application of the present formalism to a rotating magnetic dipole, the well-known 
Newtonian expressions can find general relativistic corrections that have simple physical interpretations (this will be discussed in Section|6j. 

We therefore express the electromagnetic fields in the wave-zone through two scalar functions U and V (also sometimes referred to as 
the "Debye potentials"), whose angular dependence is expanded in a series of spherical harmonics and that have a harmonic time dependence 



U(t, r, 6, 4>) = rU em (r)Y em (8, ^e" 1 "* , V(t, r, 6, 0) = rV lm {r)Y lm {9 , 0)e"^ . (86) 

Note that in the curved spacetime exterior to the star, the observed frequency for the stellar oscillations will be different for different observers 
and therefore is a function of position. In particular, if k a is the null wave-vector associated with the electromagnetic fields in the wave-zone 
and u a the 4-velocity of an observer, then to = —k a u a will be the frequency measured by such an observer. Denoting with lu r = uj{r = R) 
the angular frequency of oscillation measured by an observer at the stellar surface, the corresponding electromagnetic wave frequency at a 
generic radial position r > R will be subject to the standard gravitational redshift and given by (see Appendix |c| for a somewhat different 
derivation) 



uj(r) 



■■ N 



R~2M 
r-2M 



r 
R 



(87) 



In the asymptotically flat regions of the spacetime where we are considering the wave-zone solutions to the Maxwell equations, the electro- 
magnetic waves will have reached their asymptotic form and their frequency can then be considered to be simply u> — u> R N a = const. 

Under these assumptions, and following the approach suggested by Bouwkamp & Casimir (1954), and Casimir (1960), the components 
of the electromagnetic fields in the wave-zone assume the generic form 

1 



B r 



-Ui m Yi m e 



(88) 



B v = 



B" 



and 



r£(£+l) 

1 

r£ (£ + 1) 



d r {rUtm) doYim ^-^rVi m dsYin 

sm9 



d r {rUtm) -^—^d^Ytm + iu}rVi> m deY hl 
sins' 



tt\t~ i r \r —iujt 

& = -Vt m itm,e 



(89) 
(90) 

(91) 



E" = 



rl {I + 1) 



d r (rVim) dgY lm + ^—-rUimdcpYim 



(92) 



E' p 



dr {rVim) ^—^d^Ytm - iujrUimdeYer, 
smO 



(93) 



r£(£ + l) 

It is useful to note that while no difference would be introduced in the expressions of the radial components of the magnetic and electric 
fields [i.e. equations 1881 and <91H , the general-relativistic expressions for the remaining components would be corrected by coefficients 
of th e type ~ N or ~ 1/N, a nd effectively very small as soon as one considers regions of the spacetime away from the stellar sur- 
face JRezzolla & Ahmedovkobll) . 

In the typical configuration which we intend to consider, that is inspired by astronomical observations, the star is endowed with a 
background magnetic field that is essentially stationary on the timescale of the stellar oscillations. The background electric fields will then 
either be zero or induced by the time- variations of the background magnetic field (such as those produced by the stellar oscillations or by its 
rotation). In this case, no background electromagnetic fields will be present in the wave-zone, and expressions I88t — I93i will then effectively 
refer to the perturbations. As such they should be denoted with a symbol "5" but, in order to keep the expressions compact and since it is not 
necessary to distinguish them from the background expressions, the symbol "5" will be omitted hereafter. 
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Substituting the above expressions <88> - <93t into the Maxwell equations evaluated in a flat spacetime, these can be recast into ordinary 
wave-like equations for the Debye potentials U (or V), i.e. 

d\U -d 2 r U- -^r-Ede (sin 9d e U) - - 1 2 &U = , (94) 

H sin 6 r 2 sm 9 

which would instead take the form 

d 2 U - N 2 d r (N 2 d r u) - N 2 -^r-,dg Umedglj) + - x 2 Au 
\ I |_Hsm# \ I r 2 sin 6 

on a Schwarzschild background. Note that while equations 1941 and 1951 are not wave equations for the potentials U, V (or the corresponding 
relativistic ones U, V) , they are so for the potentials Ue m , Vim (or the corresponding relativistic ones Ue m or Vim) ■ 

The main advantage of the use of a flat background spacetime for the study of the electromagnetic fields in the wave-zone region comes 
from the fact that while no analytic solution can be found for equation I95t . an analytic solution to equation 1941 (regular everywhere except 
at r — 0), can be expressed in terms of outgoing spherical Hankel functions He(uir) in the form 

U em {r) = [£{e+l)] 1/2 H e {cjr)ui m , V lm {r) = - [£ {I + 1)] 1/2 H t (ur)vtm , (96) 

so that the components for the magnetic fields are given by the general expressions 



, (95) 



-iiot 



B r = - ^ t + 1 ^ H t (ur)ut m Yi m , (97) 

r 

B e — [ DHi(ujr)ui m deYe m - uHt(ur)vt m . " ) , (98) 

B^ - i = DHAu>r)ui m . ™ - ujH t {ior)v lm deY lrn I , (99) 

^/£(£ + 1) V sin 6» J 

while the electric field components are expressed as 



E f = -- viii±il H t {ur)vt m Y im , (100) 

r 

E e = r== DHi(ur)vt m daYi m + uH t (ur)ut m ^-^ , (101) 

E* = —i = I DHAu)r)vim . T + uH l {iur)u lm deY lm ) . (102) 

Here, we have denoted the radial derivative as DHi(cjr) = r~ 1 d r [rHi(u>r)] and we recall that the spherical Hankel functions have a simple 
radial fall-off in the case of small arguments, i.e. 

Ht(wr)ta -i(2£- l)!!^)"*" 1 , DH e (u;r) « i£(2l - l)!!(o;r)~ £ ~ 2 Lj = -I— , for cur « cu R R « 1 , (103) 

while they exhibit a typical oscillatory behaviour (in space) in the limit of large arguments (see, for instance, Arfken & Weber 2001), i.e. 

HAur) « (-i) e+1 , DHAujt) « {-if =iu>Ht, for lot -> oo . (104) 

ujt r 

We also note that the Newtonian expressions J971 - J1021 are not new and have indeed been discussed by a number of authors before 
(see, for instance, McDermott et al. 1984). However, new and general relativistic corrections are introduced in equations 1971 — J1Q21 once the 
integration coefficients ue m and vi m are specified through the matching of the electromagnetic fields 1971 — J1021 at the stellar surface with 
the help of the boundary conditions 1401 — 14 1 1 



y/I(t + 7) DHt{u) 



Ulm = 



^- / *>{*!£, [*£B* - 64K] + ^ [KK - Sv^Bi] } , (105) 

J dn {ideYt m [5v R B H - 5<S;] + ^ [8viB f R - 8v R B R ] } , (106) 



y/e(e + i)Ht(w R R)N R w B 

and of the relation 

dSl (deY lm doYl m - m Y m rnYt \ = + 

\ smfl sins J 

Indeed, and as it will become apparent in the discussion made in the following Section, all of the general relativistic effects are embodied 
in the proper specification of the integration constants 11051 — fT06l . In the case of the electromagnetic radiation produced by a rotating dipole, 
for instance, they will be responsible for the most important quantitative corrections of general relativistic nature. 
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6 A USEFUL APPLICATION: A ROTATING MAGNETIC DIPOLE 

In the following, we will discuss the application of the expressions derived in Sections [4] and [5] for a configuration frequently discussed 
in astrophysics: the electromagnetic fields and the corresponding energy losses produced when a relativistic and massive magnetic dipole 
moment rotates in vacuum. Although idealized, this simple model has been widely used to provide a first quantitative description of the 
phenomenology associated with pulsars. The expressions for the electromagnetic fields and for the energy losses through dipolar electromag- 
netic radiation have been derived long ago within a Newtonian description (Pacini 1968) and, to the best of our knowledge, have not been 
extended before to a fully general relativistic framework. This complex problem, however, can find a simple solution within the formalism 
discussed so far if it is recast into the problem of determining the electromagnetic fields in the near and wave-zone of a relativistic and 
perfectly conducting stellar crust moving with a velocity field 

1 



ou 



and whose corresponding components of the three-velocity are 



5v l 



N 



1,0,0,0 



0, 0,nrsin( 



(108) 



(109) 



where r < R and fi is the angular velocity of the star as measured by a distant observer. 

Being axisymmetric, the velocity field 11081 will not introduce first-order perturbations in the near and wave-zones if the magnetic 
dipole is aligned with the stellar rotation axis. However, if a nonzero inclination angle \ ls present between the stellar magnetic dipole and 
the rotation axis, this velocity field will introduce a time modulation of the magnetic and electric field components in the near-zone as well as 
the emission of electromagnetic waves in the wave-zone. These variations can be compared with the expressions for the electromagnetic fields 
of a rotating relativistic magnetized star available in the literature JRezzolla et alJl20 01a b; Konno & Koiima 2000), thus offering a useful 
testbed for the formalism developed so far. Before doing that, however, it is necessary to calculate the expressions for the static background 
dipolar magnetic field using the general expressions J67M69L 

6.1 Near-Zone Background Magnetic Fields of a Static Dipole 

Hereafter, we will calculate the dipolar properties of the background magnetic field by considering the general expressions <67M69> with 
I = 1, but will also allow for nonzero values of the index m to investigate the important departures from axisymmetry. Furthermore, we 
will assume the magnetic dipole to be varying on a timescale much longer than that set by the stellar oscillations so that it can effectively be 
considered static. In this case, the scalar functions Sim in 1661 are time-independent and can be easily calculated to be 



Si m (r) 



2M 3 



, *r2 2 A/ 

In N H 



after using the expression for the first Legendre function of the second kind' 5 

' x + 1 



Qi 



■In 



1 + 



1 . 



M 



Sir, 



2 \X-1 J 

Substituting J110> in equations <67t - <69| yields the following components for the background magnetic field in the near-zone 



(110) 



(HI) 



B r = 



B v = 



B" 



M 3 

N 
M 2 r 

N 



, n 2M / M 
In N 2 H 1H 



-J-lniV 2 + -i- + 1 
M N 2 



Sim' lm — 



SlmdoYl m — 



M 2 r sin6 



— lnA^ + — + 1 



M 



N 2 



Sim 0(p I 1 rn 




, „ r2 2M ( M 
lniV 2 + |H 



-£-lniV 2 + -iT + l 
M N 2 



\V2s 



sio cos6< — SKe (siie^ ) sin6H , 



[V2. 



sio sin 9 + D\e ( sue ) cos 



1 



— lnN z + —- + 1 



M 



N 2 



9k isne 



i<§t 



(112) 

J 

(113) 
(114) 



where we have chosen to consider si m = for m < and where *Re(A) is the real part of A. 

Expressions II 121 — fTT4l are the general relativistic vacuum magnetic field components in the near-zone of a relativistic spherical star 
with a dipolar magnetic field. To be fully determined, they further need the specification of the coefficients si m which can be calculated 
through the use of suitable boundary conditions across the stellar surface. As mentioned in the previous Section, these integration constants 
represent the way in which the information about the properties of the star, and their general relativistic corrections, are impressed upon the 
electromagnetic fields both in the near-zone and in the wave-zone. One way of determining these boundary conditions could be by ensuring 
the continuity of the radial component of the magnetic field across the stellar surface, but this would require knowledge of the value of the 



3 Note that in Rezzolla et al. 2001a, the incorrect expression was given for Q\, In particular, equation (89) of that paper shows the expression for Qq; the 
correct expression for Q\ was however used in the calculations reported there. 
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magnetic field components at the stellar surface and these are not necessarily easy to measure. Alternatively, the coefficients can be computed 
after multiplying both sides of J112> by Y/, m , and exploiting the orthogonality of the spherical harmonics: i.e. J dQYi m Y/, m , — 5 U i 8 mm i . 
Doing so yields the following expressions for the unknown coefficients 

M 3 r 2 



Sir. 



B r Y? m <m . 



(115) 



r 2 \nN 2 +2M(r + M) 

which must hold at any radius r and, in particular, in the asymptotically flat portion of the spacetime, where M/r <C 1 and the magnetic 
field components assume their Newtonian expressions 



(in 



2 (cos x cos 9 + sin \ sin 9 cos (f>) , cos x sin # — sin x cos # cos 0, sin \ sine 



(116) 



where, Bo = 2/i/R 3 is the (Newtonian) value of the magnetic field at the polar axis. Exploiting this property, the integrals on the right-hand- 
side of equation Jl 151 can be computed analytically giving 



Sio = 



3?r d d3 
BoR cosx 



sn = 



3tt B B 



■smx 



In the limit of flat spacetime, i.e. for M/r — ♦ and M/R — » 0, expressions il 121 — JTT4I give 



lim B r 

M/r->0,Af/«->0 



3r 3 



SljTi J 1] 



R 



b: 



(117) 



(118) 



lim 

Af/r->0,M/K->0 



4 
3T 5 



B = ^simC^Yim =--(—) d B B r R , 



1 /i? 

2 \ r 



(119) 



lim 

Af/r-»0,M/«->0 



B ,! 



3r 3 sin # 



SlmOffi I 1m 



1 /# 

2 I 7 



j^-Bfl . 



(120) 



where S^ = So (cos x cos 9 + sin x sin 6 cos 0) is the radial component of the magnetic field at the stellar surface. As expected, expressions 
ill 8t — fT20l coincide with the near-zone solutions of Muslimov & Tsygan (1986) for the dipolar magnetic field of a Newtonian magnetized 
star [cf. eqs. (18) of Muslimov & Tsygan, 1986] and that a toroidal magnetic field can be produced by a poloidal one when this is not 
axisymmetric [cf. equations II 141 and <120H . 



6.2 Near-zone Electromagnetic Fields from a Rotating Magnetic Dipole 

Given the background dipolar magnetic fields discussed in the previous Section we now examine the electromagnetic fields that are produced 
when the stellar dipole is inclined at an angle x with respect to the stellar polar axis and rotates uniformly with an angular velocity fl as 
measured by a distant observer. In this case, the integration constants st m for the background magnetic field are the same as in <117> . while 
the values for the integration constants 8se m relative to the perturbed magnetic field are computed to be (see AppendixlBlfor details) 



3SlO 



0. 



5$u 



3tt B .R 



l)sinX: 



(121) 



2 2 v 

where, again, we have chosen 8si m = for m < 0. Using now II 171 and <121> in expressions 1671 — I69i . the complete magnetic field 
components produced in the near-zone by the relativistic rotating magnetic dipole are then given by the general expressions 
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3ir 
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, *r2 2M 

In N H 
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M 



Bo (cos x cos 9 + sin x sin 6 cos A) 



(122) 



B tt = 



8M 2 r 



M N 2 



Bq (cos x sin 9 — sin x cos 6 cos A) , 



(123) 



B" 



WTN_ 

8M 2 r 



M 



N 2 



+ 1 



Bo (sin x sin A) , 



(124) 



where A = (f>—flt. Although derived in a different way, expressions 1 1 22i — fT24l coincide with the magnetic field components of a magnetized 
relativistic star in the slow-rotation approximation [cf. equations (97)-(99) of Rezzolla et al. 2001a]. Furthermore, in the case in which x = 0, 
expressions I122J — fT24l simply reduce to the magnetic field components of a relativistic dipole (Ginzburg & Ozernoy, 1964; Anderson & 
Cohen, 1970). 

We next consider the expressions for the electric field using equations 1781 — 180^ with £ = 2 for a quadrupolar electric field and exploiting 
the relation between the Xi m and Si m functions [cf. equation <77H which leaves as the only nonzero component 



X lm - y 



, *r2 2M 

IniV H 



1 + M. 



Sxim(t) 



(125) 
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With some straightforward algebra, we finally obtain that the components of the electric field produced in the near-zone by the rotating 
magnetic dipole assume the form 



3g R N R 

2 9 a N h 

j,_ 1E N 

5E -^T r n~ r 
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r \ , o 2M 2 
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sin x cos 8 sin A , 



where we express with g R the constant coefficient 

2M 2 



'''•^'-^j 111 ^-^^ 



-" -£)--" 



2 /tfV -R 

3 I i? J R-2M 



and where, in analogy with the corresponding Newtonian expression, we define the electric field Eq as 

Eo=(f H Bo)^R=(f R B )n H R. 



(126) 
(127) 
(128) 

(129) 
(130) 



It is worth noting that the coefficient f R introduced in expression <83l has also an important physical meaning since it gives a simple 
measure of the relativistic corrections which lead to an effective amplification of the magnetic fields in the vicinity of the stellar surface. In 
particular, if we indicate with Bq the radial component of the magnetic field measured at the surface of the relativistic star for an aligned 
nonrotating magnetic dipole, i.e. Bo = B r (r — R, 6 = 0, x — 0, fi = 0), then the ratio 



Bo 

B t) 



BqR a 
2fi 



= /h 



(131) 



represents the general relativistic amplification of the magnetic field strength [cf. 1122H . This amplification factor f R is a direct consequence 
of the spacetime curvature produced by the star and is shown in the small inset of Fig.Q where it takes values between 1.1 and 1.5 for the 
values of the stellar compactness MjR usually associated with compact stars. With simple algebra and proper order accounting, it is not 
difficult to show that f R — 1 in the Newtonian limit [cf. equation (93) of Rezzolla et al. 1 2001a)], so that Eo = £IRBq will then represent 
the well-known electric field induced by the rotation of a Newtonian magnetized star. 

Also for the electric fields, it is possible to show that expressions <126t - iT28l coincide with the solutions found for the exterior electric 
field of a rotating misaligned dipole in a Schwarzschild spacetime [cf. equations ( 1 24)— ( 1 26) of Rezzolla et al. 2001a]. Furthermore, in the 
limiting case of a flat spacetime, i.e. M/r, MjR — > 0, the expressions coincide with the near-zone solutions of Deutsch (1955). Overall, the 
results presented in this Section confirm the consistency of the procedure followed and show the flexibility of the formalism presented here. 



6.3 Wave-zone Electromagnetic Fields from a Rotating Magnetic Dipole 

The expressions for the electromagnetic fields produced in the wave-zone by a rotating magnetic dipole were calculated long ago in Newto- 
nian gravity, dating back to the work of Deutsch (1955), who presented them in terms of a complex but complete multipolar expansion. The 
extension of Deutsch's results to a general relativistic framework using the procedure discussed in Section[5]is rather straightforward, but we 
will consider it here only for the lowest-order multipoles. In particular, we will discuss the dipolar (i.e. those with I — 1) and quadrupolar 
(i.e. those with I = 2) parts of the electric and magnetic fields in the wave-zone since these are the ones with the slowest radial fall-off and 
are therefore the ones mainly responsible for the energy losses. 

After some lengthy but straightforward algebra in which conditions 11051 and 11061 are written explicitly in terms of the magnetic field 
components <122t - fT24l and of the velocity field J108I . it is possible to determine the different integration constants U£ m and Vim to be used 
in expressions I97i — 11021 . More specifically, when I — 1 one finds that V\ m = and the only nonzero coefficient is given by un and has an 
explicit expression 



Anx 2 3 
Mil = -JVy & H R f R Bosinx 



(132) 



Using this result, the dipolar parts of the electromagnetic fields <97> - <102l induced in the wave-zone by the rotation of the dipole are then 
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expressed as the real parts of the following solutions 

B , = _ i ^^ /sBo8in< , sinxe i[ n (r-t) + *] ( (133) 

iv R r 
B e = }_^Lf R B cose sin xe^-^ , (134) 

B i> = i$t& fBoan jMr-tn*] (135) 

2 r 

E f = , (136) 



§ i Q. R R . i[£l(r-t)+d>] ni ,„-, 

E = a Eosmve y ' vi — ±r , (137) 

2 r 

IT10 1 fi^-R . i[n(r-t)+A] D e /-no 

iv = - EocosSsmve = —B . (138) 

2 r 

Since the wave-zone is located well outside the light cylinder, i.e. atr > n c = I/O., expressions <133> - fT38l show that in this region the 
electromagnetic fields behave essentially as radially outgoing waves for which \B r /B B \ ~ \B r /B^] ~ 1/fir <C 1. 

Proceeding in a similar manner when 1 = 2 and requiring the solutions to be regular at 9 — 0, one finds that «2m = and the only 
nonzero coefficient is then given by Vaii whose explicit form is 

«2i = ^Y^^i? 5 / H B sin X ■ (139) 

As a result, the quadrupolar parts of the electromagnetic fields induced in the wave-zone by the rotation of the magnetic dipole are expressed 
as the real parts of the following solutions 

B f ' = , (140) 

i o 4 P 5 
B =---;— ~ f R B cosOsmxe , (141) 

B* = ±^-f R B (sin 2 9 - cos 2 9) sin X e l[n(r - t)+< « , (142) 

12 r 

E f = -\ ^\e sin 9 cos 9 sin xe i[n(r - t)+ * ] , (143) 

2 Jy R r 2 

E§ = iti^Eo tain 2 9 - cos 2 0) sin xe i[n(r -* )+ * ] = B+ , (144) 

12 r v ' 

1 O 3 R 4 

iv = — - — - — Eo cos 9 sin %e l J = — -B . (145) 

12 r 

Also in this case the electromagnetic fields behave essentially as radially outgoing waves for which \E r jE | ~ \E r /E^\ ~ 1/Clr <C 1. 

A few aspects of expressions 11331 — JT451 are worth noting. The first and most obvious one is that all of these components are identically 
zero when the magnetic dipole is aligned with the rotation axis (i.e. \ — 0) underlining that an axisymmetric configuration cannot emit 
electromagnetic waves. The second one is that all components have a phase term oc exp{i[fi(r — t) + </)]} expressing the wave nature of 
the solutions, where the term oc exp(ifir) is inherited from the spherical Hankel functions once u — Q. In the limit of flat spacetime, these 
solutions coincide with those in equations (87) and (88) of Michel & Li (1999) for the dipolar and quadrupolar fields, respectively. Finally, 
because expressions J133I — IT381 and 11401 — JT451 represent the dominant electromagnetic fields in the wave-zone, they provide the largest 
contributions to the electromagnetic losses which will be computed in the following Section. 

6.4 Electromagnetic Radiation Losses from a Rotating Magnetic Dipole 

Using the results of the previous Section, we can estimate the luminosity carried away by the dipolar electromagnetic radiation L em in terms 
of the integral of the radial component of the Poynting vector P 

P f dS = -!- / (exbY dS , (146) 

as 4 ?r J as V / 



LuDH e (Lur)H e (cur) [\u lm \ 2 + \v em \ 2 ] (daY tm d<,tf m + mY ^ mY ^ \ 

\ sin9 sin8 J 



(147) 
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where the integral is made over the 2-sphere <9E of radius r 3> l/f2 > R and surface element dS. 

Substituting in 11461 the expressions for the electric and magnetic fields in the wave-zone <98M99| and <10U - fT02l . it is not difficult 
to find that the radial component of the Poynting vector 

P f = _L \e $ B S ~ E 8 B $ \ = -^ r 

4tt I I M{t + 1) 

and is basically determined by the square of the integration constants ui m and vt m [cf. equations 11051 . <106H . In the case of a purely 

dipolar radiation, uii is the only nonzero coefficient [cf. equations <132H and the power 11461 radiated in the form of dipolar electromagnetic 

radiation is then 

1 1 2 o 4 C6 52 

Lern = — (\u tm \ + \Vl m \ ) = — = — -= Sin \ ■ (148) 

57T 57T DC 3 

When compared with the equivalent Newtonian expression for the rate of electromagnetic energy loss through dipolar radiation I Landa u & Lifshitzl 
ll987tlPacinll968t) . 

1 r \ Q, R B Q . 2 ,, An , 

(iemjNowt. = rr~S Sill \ , (149) 

be 3 

it is easy to realize that the general relativistic corrections emerging in expression <148> are due partly to the magnetic field amplification at 

the stellar surface [i.e. Bq — f R Bo; cf. equation I131J 1 and partly to the increase in the effective rotational angular velocity produced by the 

gravitational redshift [i.e. fl = Q R N R ; cf. equation {57}]. Overall, therefore, the presence of a curved spacetime has the effect of increasing 

the rate of energy loss through dipolar electromagnetic radiation by an amount which can be easily estimated to be 

\ 2 
jR ' , (150) 



(iem) Ncwt \N R 

and whose dependence is shown in Fig. Q with a solid line. Although expression J148> has been obtained neglecting the curvature effects 
for the electromagnetic fields in the wave-zone region (cf. Section|5j, it is straightforward to realize that already in this simplified setup, the 
Newtonian expression J149I underestimates the electromagnetic radiation losses of a factor which is between 2 and 6 in the typical range for 
the compactness of relativistic stars. Furthermore, when referred to a typical neutron star with magnetic field Bo = 10 12 G, R — 12 km, and 
angular velocity Q = 500 rad s~ , expression 11481 gives an energy-loss rate 



©• 



where f R /N R ~ 4 for M = 1.4M< 

The expression for the energy loss 11481 can also be used to determine the spin-evolution of a pulsar that is converting its rotational 
energy into electromagnetic radiation. Following the simple arguments proposed more than thirty years ago (Pacini 1968; Gunn & Ostriker 
1969), it is possible to relate the electromagnetic energy loss L em directly to the loss of rotational kinetic energy E ro t defined as 

E rot = i J d 3 x^ e -* M p(<5^) 2 , (152) 

where p is the stellar energy density. Introducing now the general relativistic moment of the inertia of the star 

1= /d 3 x^e - * (r V 2 sin 2 6), (153) 



whose Newtonian limit gives the well-known expression for a spherical distribution of matter / = (-T)Newt. = tMR 2 , the energy budget is 
then readily written as 

E Iot = j t (^inA = -L em . (154) 

Of course, in enforcing the balance 11541 we are implicitly assuming all the other losses of energy (e.g. those to gravitational waves) to be 
negligible. This can be a reasonable approximation except during the initial stages of the pulsar's life, during which the energy losses due to 
emission of gravitational radiation will dominate because of the steeper dependence on the angular velocity (i.e. -E GW oc Q 6 ). 

Expression J154I can also be written in a more useful form in terms of the pulsar's most important observables: the period P and its 
time derivative P = dP/dt. In this case, in fact, using expression 11481 and J154I . it is not difficult to show that 



*\ i 

iteJl969h 



where the Newtonian expression is given by ( Gunn & Ostrike 

(P*W ^ (g) *f . (156) 

Also in this case it is not difficult to realize that general relativistic corrections will be introduced through the amplification of the magnetic 
field and of the stellar angular velocity, as well as of the stellar moment of inertia (i.e. I /I). 

Expressions 11481 and 11551 could be used to investigate the rotational evolution of magnetized neutron stars with predominant dipolar 
magnetic field anchored in the crust, either when isolated JPage et al. 2000) or when present in a binary system JLavagetto et a l. 2004). Indeed, 
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Figure 1. General relativistic amplification of the energy loss through dipolar electromagnetic radiation as a function of stellar compactness. The solid line 
refers to the result derived here [cf. eqs. 11481 and 11491 1. while the dashed line refers to the phenomenological estimate made by Page et al. (2000) (PGZ). 
Shown in the inset is the magnetic field amplification produced by the spacetime curvature for different values of stellar compactness [cf. eq. 11311 1. 



a detailed and accurate investigation of this type has already been performed by Page et al. (2000), who have paid special attention to the 
general relativistic corrections that need to be included for a correct modelling of the thermal evolution but also of the magnetic and rotational 
evolution. It should remarked, however, that in their treatment Page et al. (2000) have adopted an expression for the loss of rotational kinetic 
energy which is similar to expression 11551 in that it accounts for the magnetic field amplification due to the curved background spacetime, 
but that does not include the corrections due to the gravitational redshift. As a result, the general relativistic electromagnetic luminosity 
estimated by Page et al. (2000) is smaller than the one computed here and it is shown with a dashed line in Fig.Qfor comparison (labelled as 
PGZ). As the authors themselves discuss, their choice was based on simple phenomenological considerations and was made in the absence 
of a more systematic treatment of the electromagnetic fields produced by a rotating and magnetized relativistic star. We hope that the results 
reported here will be of help in all of those investigations that consider the thermal, magnetic and rotational evolution of magnetized neutron 
stars. 



7 CONCLUSIONS 

There is undoubted observational evidence that strong electromagnetic fields are often present in the vicinity of relativistic stars. At the same 
time, it is natural to expect that these stars will be subject to perturbations of various types which will induce oscillations responsible for 
the emission of both electromagnetic and gravitational waves. Clearly, the possibility of detecting simultaneously both types of signal will 
provide important information about the mass and radius of these objects and hence give the ability to investigate the properties of matter 
at nuclear densities. Building on the general relativistic treatment of electrodynamics of magnetized neutron stars developed by Rezzolla 
et al. (2001a,b), we have studied here the electromagnetic fields generated when the perfectly conducting crust of a magnetized relativistic 
star possesses a nonzero velocity field. The star has been modelled as a relativistic polytrope with infinite conductivity and in vacuum. The 
relativistic corrections induced by a rotation in the spacetime as well as those produced by the electromagnetic energy densities have been 
neglected. We expect both of these app roximations to be satisfa ctory for the large majority of rapidly rotating and magnetized neutron stars. 
As in our previous investigations JRezzolla et all2001albh . we have paid special care to providing expressions for the electromagnetic 
fields that besides being analytical, are written in a form which is as simple as possible, that highlights the relation between the general 
relativistic and Newtonian expressions, and that can therefore be used also in astrophysical applications. This first paper has been mostly 
devoted to the development of the mathematical formalism which distinguishes the solution near the stellar surface (where the electromagnetic 
fields are quasi-stationary), from the one at very large distances (where the electromagnetic fields approach a wave-like behaviour). In both 
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of these regions and with suitable approximations, in fact, exact solutions to the general relativistic Maxwell equations can be found and 
have shown to contain two main corrections when compared with the corresponding Newtonian expressions. The first one is related to the 
amplification of the electromagnetic fields on the stellar surface caused by the background spacetime curvature, while the second one is due 
to the gravitational redshift which will affect the electromagnetic waves as these propagate in the curved spacetime. It should be noted, in 
fact, that while the wave-zone solutions to the Maxwell equations formally coincide with the ones derived in a flat spacetime, integration 
constants defined at the surface of the star, imprint important general relativistic corrections which, in turn, convey information about the 
mass and radius of the oscillating star. 

The general expressions for the electromagnetic fields presented here do not refer to any specific magnetic field topology or velocity 
field and are, in this sense, completely general. However, to validate the expressions derived and to consider an astrophysically important 
testbed, we have also examined the form which these electromagnetic fields assume when the background magnetic fields are of dipolar 
type and the crustal velocity field is that produced by the uniform rotation of the dipole. This configuration is often used as the simplest 
model for the electromagnetic emission observed from pulsars and has been extensively investigated in the past. While analytic expressions 
are available within a Newtonian framework both in the near-zone and in the wave-zone, the general relativistic expressions are limited to 
the near-zone. A direct comparison with these expressions and the perfect agreement found have shown that the formalism developed here 
is flexible and robust. Furthermore, the use of analytic expressions for the components of the electromagnetic fields in the wave-zone has 
allowed to calculate, for the first time, the general relativistic expression for the electromagnetic energy loss through dipolar radiation. A 
simple comparison with the corresponding Newtonian expression shows that the latter under-estimates this loss by a factor of 2-6, depending 
on the stellar compactness. 

Work is now in progress to apply the formalism presented here to those spheroidal and toroidal velocity fields most frequently encoun- 
tered in the study of stellar oscillations. The results will be presented in a forthcoming companion paper JRezzolla & A hmedov 2004). 
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APPENDIX A: DERIVATIVES OF LEGENDRE FUNCTIONS 



In this Appendix we briefly sketch the lenghty but simple algebra needed for the calculation of the functions Qe m and their radial derivatives, 
and which are used in the calculation of the electric field in the near-zone. Using the recurrence formulae (see, for instance, Jeffrey 1995) 



' dx 



(Al) 



(I + 1) Qi+i - (21 + 1) xQi + lQ e -! = 



(A2) 



where, we recall, x = 1 — r/M, it is possible to show that the radial second order derivative of Legendre functions of the second kind takes 
the form 
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As a result, the expressions contained in the electric field components 1781 — 1801 can be written explicitly as 
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APPENDIX B: CALCULATION OF THE INTEGRATION CONSTANTS d T 5Si M 
Inserting into 1701 the value of the magnetic field 1471 - 1491 at the surface of the star 

B% = -±(VlS)\ r=R , B° = 
and making use of expression I55i . it is possible to derive that 
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(B2) 



where V n is the angular part of the Laplacian, i.e. 



1 1 2 

-^—rde (sin^fl) + . %- di 
sm 6 sm 6 



(B3) 



It is interesting to note that in the limit of flat spacetime our expression <B2t coincides with equation (B2) of Timokhin et al. 2000, 
which was instead derived after requiring the continuity of the tangential electric field at the stellar surface. Multiplying this equation by 
Yi m , integrating it over the solid angle and using equation I66i . we can finally obtain an expression for the coefficients dt{SSe m ) in the 
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expansion of SS in spherical harmonics 

dt(SSim)\r=R = jj {- (l + jjl\ Ql + Ql-1 } <■>,**■„. (t) 
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APPENDIX C: A DIFFERENT APPROACH TO THE GRAVITATIONAL REDSHIFT 

We here provide a somewhat different derivation of the well-known result that photons in a gravitational potential undergo a gravitational 
redshift. For doing this we will exploit the existence in this spacetime of a timelike Killing vector £ Q such that 

£a;/3 +&;<» = . (CI) 

Electromagnetic waves propagate along geodesies and the associated null wave-vector k a will be tangent to these trajectories and 
parallel transported along them, i.e. 

k a .^ = . (C2) 

We now introduce two important frequencies: the first one is the frequency measured by an observer with four-velocity u a and defined as 

u) = —k a u a , (C3) 

while the second one is the frequency associated with the timelike Killing vector £ a and defined as 

cut = -k a £ a . (C4) 

The two frequencies have marked differences: while IC3i depends on the observer chosen and is therefore a function of position, <C4> 
is a conserved quantity that remains constant along the trajectory followed by the electromagnetic wave. This is easily verified after using 
equations <CU and <C2> which give 

(fc Q C Q ), fl k = (r^) ;)3 k 1 * = k a . kPt a + U^k" = . (C5) 

We can now exploit this property to measure how the frequency changes with position and is therefore redshifted in the spacetime J2j. 
Assume the Killing vector to have components 

C =( 1,0,0,0 ) ; ksW -1,0,0,0], (C6) 



so that lo^ = — feo =const. The frequency of an electromagnetic wave emitted at the surface of the star r = R and measured by an observer 
with four-velocity u a parallel to £ a (i.e. a static observer) will be 

Lu R =-(k a u a )\ r = R = N- 1 Lu i , (CI) 

so that at a generic radial position r 

N(r)uj(r) = const. = N R u R . (C8) 

This expression coincides with equation 1871 in the text. 

APPENDIX D: MAXWELL EQUATIONS IN A SPACETIME WITH A ROTATION 

If the spacetime admits a rotational Killing three-vector j3, the full set of Maxwell equations can then be written as (Thorne et al. 1986, 
Rezzolla et al. 2001a, b) 

V • B = , (Dl) 

V ■ E = 47rp e , (D2) 
(dt - C )B = -Vx (NE) , (D3) 
(d t -£ g )E: = -Vx(NB)-4irJ, (D4) 
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where, again, V represents the covariant derivative with respect to spatial part of the metric, C s is the Lie derivative along /?, and N is the 
lapse. The system of equations JDU - JD41 is completed with Ohm's law, which takes the general form 

j =a\E+ (v + 0) x fi] . (D5) 

The Lie derivative of a three-vector A along the vector field j3 is not particularly difficult to calculate. However, the form of its com- 
ponents in a spacetime with line element J2j and a frame {6j-{5J is not easy to find in the literature. For this reason we give them here as a 
useful reference 
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If the Lie derivative is now taken along the rotational Killing vector 

0= \-^ L - a e<t, = -ursmd e , (D9) 

2 r sin 

as for the spacetime of a rapidly rotating relativistic star, expressions <D6MD8l reduce to 

[CpAX = -uB^A* , (D10) 



CpA) =-Lod<t,A e , (Dll) 
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